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ABSTRACT 



We study the large N limit of the Itzykson - Zuber integral and show that 
the leading term is given by the exponent of an action functional for the complex 
inviscid Burgers (Hopf) equation evaluated on its particular classical solution; the 
eigenvalue densities that enter in the IZ integral being the imaginary parts of 
the boundary values of this solution. We show how this result can be applied to 
"induced QCD" with an arbitrary potential U(x). We find that for a nonsingular 
U (x) in one dimension the eigenvalue density p(x) at the saddle point is the solution 
of the functional equation G+(G-(x)) = G-(G+(x)) = x, where G±(x) — ^U'(x)± 
mp{x). As an illustration we present a number of new particular solutions of the 
c = 1 matrix model on a discrete real line. 



1. Introduction 



The analysis of most any multi-matrix model is based on the Itzyk- 
son - Zuber formula [1] 

Hi, v) = f vUe^Wl = det[eArtol (1 1) 

where (f) and x are arbitrary N x N Hermitian matrices, A(0) denotes 
the Van der Monde determinant constructed out of the eigenvalues 4>i 
of these matrices, A(0) = n^<j(^ — and P[/ is the Haar measure 
on the unitary group U(N). The way to apply it is to diagonalize the 
matrices entering the model and to carry out the integration over the 
diagonalizing unitary matrices using (1.1). Then, as N — > oo, the resid- 
ual integrations (over the eigenvalues of the matrices) can be treated 
by the saddle point method. In this limit the eigenvalues at the saddle 
point can be described by their distribution on the real axis: if di is the 
number of eigenvalues between and (f) + d(f), then there is a smooth 
function p{4>) such that di/N = p((f))d(f). If we keep the distribution 
p{4>) of the eigenvalues 4>i entering (1.1) (as well as a(x) of Xi) fixed, 
the integral (1.1) posesses the N — > oo asymptotic expansion, 

7(0, x) ~ exp{7V 2 F o (0, X ) + FM, X ) + N- 2 F 2 (<f>, X ) + •••}• (1-2) 

Nevertheless, although there is an exact explicit formula (1.1) for the 
Itzykson - Zuber integral for any finite N, the structure of the coefficients 
F n is not known yet. For many applications it would be important to 
understand it. 

An interesting example of such an application is the so-called "in- 
duced QCD" [2] which is the lattice gauge theory of a scalar matrix field 
$>(x) defined by the action 



X 



S = Y,NTr -U(S(x)) + ^(x)U^(x)$(x + /ia)Ut(x) 



A* 



;i.3) 



(fi runs over the 2D lattice vectors on a hypercubic lattice of dimension 
D, the field $(x) transforms according to the adjoint representation of 
the gauge group U (N)). It describes the high temperature limit of D + l 
dimensional QCD and, in addition, can be related to the matrix models 
of c = D strings. The Itzykson - Zuber integral turns out to be the 
main computational tool necessary to analyze this theory. 



1 

2 



In this paper we claim that large N limit of the IZ integral is given 

by 

/(0, X )-exp(iV 2 Fo(0,x)), 

Fo{<f>,X) = S[p,cr] + J p(x)x 2 dx + J a{y)y 2 dy^ 

[ p{x\)p{x2) In \x\ — x<2\dx\dx<i (1-4) 
+ J (j(yi)a(y 2 )ln\yi -y 2 \dyidy2^ 
where S[p, a] is the classical action corresponding to the Hopf equation 



calculated on its solution satisfying the boundary conditions 

Imf(x,t = 0) =irp(x) 
lmf(x,t=l)=7rcr(x) 



We will show that the action S[p,a] defines the field theory with the 
Hamiltonian 

H[p,n] = \J ^)j(^HM) 2 - ^Ax)\dx (i.7) 



where U(x) is the canonical momentum conjugate to p(x). In this no- 
tation 

f^t) = ^^-+mp{x,t). (1.8) 

Quite remarkably, (1.7) is the nonunitary version of the Hamiltonian, 
arising in the collective field formulism of c = 1 string theory [8] . 



The plan of this paper is as follows. In sections 2 and 3 we establish 
the relation between the IZ integral and the Hopf equation, and prove 
the boundary conditions (1.6). In section 4 we apply our results to 
analyze "induced QCD" with an arbitrary nonsingular potential U{x). 
We obtain that in the one-dimensional case the eigenvalue density p(x) 
as a solution of the master field equation must satisfy the functional 
constraint G+(G-(x)) = G-(G+(x)) = x, where G±{x) = ^U'(x) ± 
mp(x). Although we were not able to find an explicit general solution 
to this equation, we show how one can generate a number of particular 
solutions of it. We use them to analyze the c = 1 matrix model on a 
lattice for finite lattice spacings. 



One of the ways to approach the problem of the large N expansion is 
to write the differential equation for the IZ integral and specify the cor- 
responding boundary conditions. So, first we will derive the differential 
equation and then identify the relevant boundary conditions. 

It is convenient to consider the one parameter family of integrals : 



so that I((f),x\t) = exp(-gTr(0 2 + x 2 ))l(^, ^) and the Itzykson - 

Zuber integral per se corresponds to t = 1. Since I((f), x\t) depends only 
on the eigenvalues of and x, one can show [1] that (2.1) satisfies the 
partial differential equation 



where A(0) = n i<<7 -(0i - (f)j). 

This equation is the linear heat equation for the quantity 7(0, x\t) = 
A (0)7(0, x\t) and can be solved exactly to obtain the Itzykson - Zuber 

formula (1.1). However, unlike 7(0, x|t), I(4>,x\t) does not have any 
regular large N asymptotics. Therefore we have to study the equation 
for I. Let us use the ansatz 



2. Reduction of the Itzykson — 
Zuber integral to the Hopf problem 




(2.1) 




(2.2) 



7(0,x|f)=exp(wV[p,<7|*]) 



so that from (2.2) follows 



aw lia 2 ,,, *,^fdw\ 2 „-^„, ls dw 

1 = 1 1 1=1 ' 1 = 1 



where 

The equations above are exact. Now we assume that W[p, a\t] has 
a smooth large N limit, i.e. it is a smooth functional of the eigenvalue 
densities p and a. This implies that ~ (9(1 /TV) and 



1 / 5 <W 



50i N\dx5p(x) 



(2.4) 



Thus we can neglect the first term on the left-hand side of (2.3) (it is 
0(1/ N) while the two other terms are 0(1) )to get 



<9iy J^/dJy\ 2 dW 



Substituting 



H/ = s - wi EH<a< - fcl - ^2 E ln l» - %l ( 2 -6) 

we find 



* If we would have been looking for W such that I = exp(N 2 W), we would not 
have been able to neelect the first term. 



Finally, we use the identity ± Y,iLi ^ 2 (0i) = w ^i^j (^-<fe) 2 t0 obtain 



9^- ATSrf dS Y 1 1 



It is very useful to rewrite this equation in terms of the eigenvalue 
distributions. To this end we replace all the sums by the integrals ac- 
cording to the rule 



1 N f 
— -> / p(x)dx. 



i=l 



However, the sum Yli^j ^.l^.y needs to be treated specially, since 

the relevant integral J J j^ P yp dy diverges. In fact the 0(1) contri- 
bution to this sum comes from the pairs of the eigenvalues <pi and <pj 
such that \i — j\ «iV. Thus fa — (f)j « \i — j\/ '(iV p(fa)) and 



Replacing sums by integrals and using (2.4) we see that the continuum 
version of (2.7) is 



2 £ = (s^)'-T*4 (2 ' 8) 

This is a nonlinear differential equation for the functional S[p(x), a(y)\t]. 
It would be extremely hard, if possible at all, to find its general solution. 
Fortunately, we have to find only a particular solution - the one which 
corresponds to the large N limit of the IZ integral. To do so we will 
guess what the desired solution is and then prove the conjecture. 



First note that (2.8) can be represented in the Hamilton - Jacobi 
form 




where r = — t, 




(2.10) 



is the Hamiltonian , and U(x) - the canonical momentum conjugate to 
p{x). Our conjecture is that S[p(x), a(y)\r] is given by the action of 
the dynamical system defined by (2.10) calculated along the trajectory 
which connects the two points p(x) and a(y) in the configuration space 
of this system. The trajectory has to satisfy the equations of motion 
and take p(x) to a(y) within the time interval r. 

By construction the so defined S[p(x), a(y)\r] solves (2.9). One still 
has to show that it gives the right large N asymptotics of the IZ integral. 
This will be done in the next section. In the meantime we will investigate 
the properties of this solution in more detail. 

To find S[(/>, x] we have to solve the Hamilton equations following 
from (2.10). These are 

dp(x, t) 
dt 

dU(x,t) 
dt 



d_(, t) du( x ,ty 

dx\ ' dx 



1 / gnQzyQ 

2 { dx 



7T 



Introducing 



we get 



* The discrete version of this Hamiltonian, the one corresponding to (2.7), de- 
scribes the completely inteerable Caloeero theory [91. 



dp 9 . N 

f a 2 » • (2-12) 

These are the Euler equations for a one dimensional liquid obeying the 

equation of state P = —^p (x). They can be easily solved. Indeed, let 
us introduce a new unknown function 



f(x, t) = v(x, t) + i7ip(x, t). 
Then / satisfies the Hopf equation 



df d 



Now we have to specify the boundary conditions for (2.12), (2.13) 
that lead uniquely to the original problem (1.1), (1.2). Our conjecture is 
to identify the endpoints of the classical solution we are looking for with 
the eigenvalue densities of the matrices (j) and x- Therefore we demand 
that 

J p(x, t = 0) = p(x) 
[ p(x, t = 1) = a(x) 

or, equivalent ly, 

J Im/(x, t = 0) = 7rp(x) , . 

\Imf(x,t=l)=ira(x) ' [Z ^ V 



Finally, taking into account that due to (2.1) Fq((J), \) = W[p, a\t] + 
2^Tr(0 2 +x 2 ) and remembering the definition of S[p, a], (2.6), we obtain 
our final result (1.4), (1.5), (1.6). 

Is it possible to find an explicit formula for Fo(0, x) ? To do this we 
have to find the trajectory obeying (2.13), (2.14). The Hopf equation is, 
in principle, exactly soluble. For a real function / the general solution 
of (1.5) can be written down in the parametrized form 



x =R(a) + F(a)t 
f(x,t) =F(a) 



(2.15) 



where R(a) and F(a) are the two functions to be determined from the 
initial conditions supplementing (2.13). In fact, because of the invari- 
ance with respect to reparametrizations of a, (2.15) depends on only 
one function, consistent with the fact that the Hopf equation is of first 
order. 

Our problem is a little more complicated, because / is complex 
valued. One can generalize the parametric solution to describe this 
situation provided that the initial condition fo(x) = f(x,t = 0) can be 
analytically continued to the whole complex plane as a function of x. 
This is indeed true in many applications. 

In this case F(a) and R(a) can be viewed as the analytic functions 
of the complex parameter a. To impose the boundary conditions we 
introduce another pair of functions, 

G + (x) = x + f(x, t = 0), G-(x) = x — f(x, t = 1) 

and use (2.15) to get 

G+(x) =x+[Fo R-^ix) = [(F + R)o R-^ix) 

as well as 

G-(x) = x-[Fo(F + i?) _1 ](x) = [R o (F + R)' 1 }^) 

so that 

G+(G-(x)) = G-(G+(x)) =x 
ImG+(x) = 7rp(x) (2.16) 
ImG-(x) = — 7ra(x) 

where o denotes functional composition. Given p{x) and a(x) we must 
solve (2.16) for, say, G+(x), then the full initial condition for the Hopf 
equations (2.13) will be given by f(x,t = 0) = G+(x) — x. This allows 
one, in principle, to calculate the trajectory connecting p(x) and a(x) 
and the action along it, which, in turn, gives the large N asymptotics of 
the Itzykson - Zuber integral. 

We were not able to obtain an explicit expression for Fq in (1.2) as 
a functional of p and a. However, in some cases (2.16) already provides 
enough information. This is the case in "induced QCD" . We will return 
to this question in section 4. 



3. Boundary Conditions 

In this section we will prove our conjecture about the boundary 
conditions for the Hopf equation. We will show that the quantity S[p, a] 
entering (1.4) is indeed the classical action along the trajectory of (1.7) 
with the endpoints 

p(x, t = 0) = p(x); p(x, t = 1) = a(x). (3.1) 

To this end we note that the Itzykson - Zuber integral 



7(0, X ) = J VUe NT ^x^] (3 . 2; 



satisfies the set of identities with higher order differential operators 



n — ll \Nd6) 



H q N I{^ X )= [Tr(x 9 )]/(0,X), k = l,...,N. (3.3) 

Since 7(0, x) depends only on the eigenvalues of the matrices (f) and X-, 
these identities reduce to 

V q I(<f> lX ) = [MX Q )]I(^X) (3.4) 

where 

1 N 8* 

To take the large N limit of the differential operator in (3.5) we 
expand it explicitly to obtain 

q\^_^(ki-l)ni+...+(k s -l)n s 



Mk 7 ? 1 . . . A#«m! . . . n s \(q - A - mh - ... - n s k s )\ 

K\ ,. . . ,K S 

X;m,...,n s 
n 

/ o \ q-(X+n\ki + ...+n s k s ) 
A/„ , , \\ni /ry , , \\n s I O \ 



(v^ n ))\z kl (0 n )y\..(z ks (0 n )y 



dcj) 



n , 



(3.6) 

where 



k^n 



On - (f>k) 



In the large N limit 



V^n) -> N V [ = O(N) 



'n 



X 



where V denotes the principal value prescription, and, naively, Z m ((f) n ) 
would go like O(N) as well (since it includes one summation over k). 
However, it is not true because of the singularity near (f) n ~ (f)^ where 
0n - <^k ~ {n - k)/ (Np((f> n )) so that 



_ ( 2N m ({m)p m {(j) n ) + 0(N m - 1 ), if m is even 
" \ 0(N m - l )^m is odd 



(3.7) 



where 



oo 1 
C(m) = V — 

^ v ' / j fern 



k=l 



is the Riemann zeta function. 

When the operators V q act on 7(0, x) = exp(iV 2 Fo[/9, cr]), the lead- 
ing TV — > oo contribution comes from 



d 



v 



d(f) 



exp(iV 2 F ) = 



n , 



n p( n ^A p + G(N p - 1 ) 



exp(7V 2 F ) 



and is given by 



V q exp(N 2 F ) = lw Y, 



ki,...,k s even 

all dif ferent,all>2 

X;ni,...,n s 
n 



q\(_^(ki-l)rn+...+(k s -l) 

X\k^...k^n 1 \...n s \ 



1 



(q - A - mki - ... - n s k s )\ 

(d 6Fq \q-(X+niki+...+n s k 
dx5p(x) x=4> n ) 



i=i 



(3.8) 



+ Oim- 1 ) } exp{N 2 F ) 



Performing the summation over hi, . . . , k s ; ni, . . . , n s and s; using 



1 N f 

i=l 



as TV — > oo we obtain 



P,exp(AT 2 Fo[p,<7]) = 
1 1 



dx 



2mq + l 
/ a 5F [p,a 
\dx 5p(x) 



dx 5p(x) 



-V{x) — l7Tp(x) 



9+1' 



exp(AT 2 F [p,a]) 



a{y)y q dy 



exp(N 2 F [p,a}) 



(3.9) 



and thus the large N limit of the identity (3.3) is given by (3.9). If we 



* To do the sum it is convenient to introduce an operator = Yl^=o ^~^ > Q- 
Then the summation can be performed using the identity 

sinh(7T2;p(0)) 



exp 

with the result 



-2 E 



nzp(<p) 



Poo exp(iV 2 Fo) = — / ttyexp 

7TZ 



now substitute ansatz (1.4) into (3.9) we will get that the quantity 

5S[p, a] 



U(x) = 



5p(x) 



should satisfy 
1 1 



2m q + l 



dx 



—Ii{x) + V{x) + mp{x) 



d 



U(x) + V{x) — iTrp(x) 



dx 

= J °{y)y q dy 



(3.10) 



To prove our conjecture we have to show that the action S[p,a] on 
the trajectory satisfying the boundary conditions (3.1) does have the 
property (3.10) as a consequence of the equations of motion (2.12). To 
see that this is indeed the case, we introduce the set of time dependent 
functions on the phase space (here v(x) = J^U(x)) 



H q (U, p, t) = 2mt( ^ ~ j-j J dx^(x- (v(x) - lTTp{x))l > j 



— (x— (y(x) + mp{x))t^ 



q+l 

q +r 



(3.11) 



Using (2.12), it is easy to check that 



dH q 
~dt 



= 



on the equations of motion of our original Hamiltonian (1.7). So we 
have found an infinite number of conservation laws corresponding to 
the problem. Moreover, 



H q {IL,p,t = l) = 



1 



- / dx< \x — (v(x) — iTip(x)) ) 
2iri(q+l)J \\ v w iK n ) 

— [x— (v(x) + inp(x))^j 



q+l 
q+l 



and 



H q (U,p,t = 0)= / a(y)y q dy. 



Now, since H q are the constants of motion, we conclude that (3.10) does 
hold, and, therefore, (1.4) is proved. 



4. Induced QCD With an Arbitrary Potential 

Induced QCD is the model defined by the partition function [2] 

Z= fv^x)VUp(x)exp^2NTr[-U(${x)) + 

J X 

J2 <S>(x)U^x)<S>(x + fm)Ul(x)] . 

A* 

The link fields U^{x) can be integrated out giving the effective theory 
of the field $>(x), with the action 



S=Y,[-NU(Ux))} +E In I'M*) "<fe(*)l 



x,i 



+y^ln 7(0(x), 4>{x + //a)) 



(4.1) 



In the large N limit the partition function of this effective theory is 
dominated by a translationally invariant saddle point which obeys the 
equation 



N ^ <f>i - ^ 



(4.2) 



<f>=x 



Using (1.4) and N§- = (£Jgj) | (^n(x)) | . = „(*), one 



dx8p(x)/ \x- 

can rewrite it in the form 



or, equivalent ly, 



»(*) = 2^C/ (x) - x + — V J — . (4.3) 

Migdal [3] has converted this into a nonlinear integral equation for p{x). 
On the other hand, from the results of section 2 it follows that the 
problem can be reduced to a differential equation. Indeed, 

v(x) = Ref(x,t = 0) = -Ref(x,t = 1) 

where f(x,t) is the solution of the boundary problem 

dt J dx J (4.4) 
Imf(x, t = 0) = Imf(x, t = 1) = irp(x) 

Equations (4.3) and (4.4) provide a description of induced QCD at large 
N. 

Let us study the simplest case D = 1. Then for a nonsingular 
potential U (x) 



G+{x) = x + /(x) = x + f (x) + iirp(x) = -[/'(x) + iirp(x), (4.5) 

and we obtain the constraints determining the saddlepoint of one di- 
mensional induced QCD in the form 

r G+(G-(x)) = G-(G+(x)) = x 

< ReG±{x) = \U'{x) . (4.6) 

[ \mG±{x) = ±7rp(x) 



Note that in this case G+(x) = G-(x). Moreover, since D = 1 induced 
QCD is equivalent to the Hermitian matrix model on a discretized line, 
one can use (4.6) to investigate the properties of this model. 



* Let us, however, notice that for some singular potentials (say, those considered 
in [4]) the analyticity hypothesis (2.15) does not hold and one has to study the 
general equations (4.3) and (4.4). 



In fact, it is easy to construct a number of exact solutions to (4.6) 
using the following device. Consider an arbitrary real analytic symmet- 
ric function P(z,w) = P(w,z) (say, a polynomial) of two variables, z 
and w. If we define and G- as the appropriate roots of 

P{x,G+(x)) = P{G-{x),x) = 0, 

then, by construction, G+(G-(x)) = x, and the roots , G- can be 
chosen to be complex conjugate, i.e. be the solutions of (4.6). To clarify 
this point, let us consider an example 

P(x, y) = x 2 + y 2 — m 2 xy + p. (4-7) 

For this P(x, y) 



m 2 / LL^X 2 

G±(x) = —x±i]Jp — 

/ 2 2 

so that 7rp(x) = y \i — ^j- (where \i is fixed by J p{x)dx = 1 to be 

p = V m 4 — 4) -the semicircular law, which is known to solve induced 
QCD with the quadratic potential U(x) = m 2 x 2 /2 [5] . One can verify 
explicitly that (4.6) indeed holds for these G±. 

To obtain a more general solution, we can perturb (4.7) by mono- 
mials of higher order. Then the relation G+(G-(x)) = G-(G+(x)) = x 
will not be violated, at least for small enough coefficients in front of 
these monomials. The simplest generalization is 

P(x, y) = x 2 + y 2 — bx 2 y 2 — m 2 xy + ^{x + y) + a. 

The corresponding potential and density are given by 

m 2 x — £ 



U\x) = 

± — ux~ 

7rp(x) = - — ~^~~2 \/ (1 — bx 2 )(x 2 + £x + a) — ^(m 2 x — ^) 2 



(4.8) 



To see the physical meaning of this solution, let us consider its contin- 
uum limit. In the continuum limit we must scale 



m 2 = 2 + a 2 M 2 ;b = Ba 3 ; 



^ = Ea 3/2 ; U($) = <S> 2 + aW(tp) , 

where a is the lattice spacing. As a — > 0, the partition function of 
induced QCD reduces to 



Z = 



J p (/? (^) e -^/^[Tr{^^W 2 +W(^(t)) + |(a^(t)) 2 }] ( 4-9 ) 



In our example 

117// \ ivi y> -\- z£>(^ _ 
^ = 1 - BaV 

We see that as a — > 0, 



VF(c^) -> — + — - Hp 



which is a general quartic potential. Note that our solution applies not 
only in the continuum limit, but for finite lattice spacings as well. 

Therefore, it is possible to find nontrivial potentials, which have 
a well-defined continuum limit and, at the same time, admit an exact 
solution of the model. 



5. Conclusions 



We have shown that the large N behaviour of the Itzykson - Zuber 
integral, which, by itself, is not dominated by any of its N\ saddlepoints 
on the group manifold U (N), is nevertheless controled by a well defined 
classical trajectory in the space of eigenvalues. We obtained the repre- 
sentation of its asymptotics in terms of that trajectory and reduced the 
problem to a set of two algebraic constraints. It turned out that the 
analysis of induced QCD is the problem, dual to the computation of the 
Itzykson - Zuber integral: to solve induced QCD we have to find p(x), 
given U (x) in (4.6), while in order to compute the integral we must solve 
for U(x), when p(x) is known. We demonstrated how some particular 
solutions of this problem can be obtained. 

Finally, it is necessary to make the following remark. Makeenko [6] 
claimed that induced QCD in D dimensions with an even potential U (x) 
has the same saddle point as the one-matrix model with the potential 

U (x) related to U (x) by 



However, it appears that there are some counterexamples to this state- 
ment. Let us, say, consider the D = 1 case. Choose 

2 2 4 

i , m x qx 

It is easy to find the solution of this model up to the first order in g. One 
can do it using (4.8) with £ = and b = ^ <C 1. Then the eigenvalue 
density p(x) reduces to 



U'{x) 



2D - 1 



1 





TT 



where 



9 



1 (m 4 -2) + 0(g 2 ), d 



c = 



2pm 



and p = Vm 4 — 4. Then the constant a is fixed by the requirement 
J p(x)dx = 1. It can be checked that the same result follows indepen- 
dently from the perturbation theory for MigdaVs equation [3]. However, 



the results following from (5.1) do not agree with that. Indeed, the 
corresponding one-matrix potential would be 

c%) = ^ + ^ 4 + o( fl 2 ). 

This matrix model is easy to solve exactly [7] . Expanding the result in 
g one gets 

p{x) = — (cx 2 + d) \ja — x 2 , 

7T 

with 

We see that c ^ c already in the first order of perturbation theory, 
indicating problems with (5.1). 
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